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Abstract 


Some  thermodynamic  relations  are  derived  along  the  principal  Hugoniot  of  materials  for 
which  the  Griineisen  relation  is  a  function  of  volume  only.  Rather  than  being  expressed  in  terms 
of  traditional  thermodynamic  variables,  such  as  volume  and  temperature,  the  relations  are 
expressed  in  terms  of  the  shock-Hugoniot  behavior  and  of  a  term  grouping  that  is  related  to  the 
Griineisen  function.  By  so  doing,  a  new  perspective  is  gained  on  both  the  nature  of  the  Hugoniot 
as  well  as  the  interrelation  of  thermodynamic  quantities  along  the  Hugoniot. 
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1.  Background 

Thermodynamic  relations  are  traditionally  expressed  in  terms  of  standard  thermodynamic 
state  variables,  such  as  volume  and  temperature.  In  this  report,  we  examine  the  thermodynamics 
of  solids  along  the  principal  Hugoniot  of  the  material.  We  break  from  the  traditional  approach, 
however,  and  express  the  behavior  of  thermodynamic  paths  through  the  Hugoniot  [in  the 
pressure-volume  (p,V)  and  energy-volume  (E,V)  planes]  in  terms  of  the  the  behavior  of  the 
Hugoniot  itself  and  a  term  grouping  that  is  related  to  the  Griineisen  function,  T  -  V(dp/dE)v.  In 
this  study,  we  consider  only  materials  for  which  T  -  T(V)  alone. 

The  equation  of  state  traditionally  expresses  pressure  in  terms  of  volume  and  temperature, 
or  alternately  in  terms  of  volume  and  internal  energy  (if  the  specific  heat  behavior  is  known). 
For  the  modeling  of  solids,  the  so-called  Griineisen  form,  whereby  T  -  T(V )  and  p  -  p(E,V),  has 
been  widely  accepted  to  model  material  behavior  (in  the  absence  of  phase  changes)  over 
pressures  into  the  megabar  range.  In  this  form,  the  pressure  and  energy  at  arbitrary  states  are 
related  back  to  conditions  along  a  reference  path  of  known  thermodynamic  characteristics 
(denoted  by  subscript  “ref'): 


P-Pref-  >  (!) 

where  the  variable  \|r  is  given  by  \|f  -  V/T  -  ( dE/dp)v .  In  this  equation,  E  and  V  denote  the 
specific  energy  and  specific  volume,  respectively.  If  the  equation-of-state  reference  curve  is 
taken  as  the  principal  Hugoniot  (the  Hugoniot  originating  at  zero  temperature  and  pressure), 
where  the  Rankine-Hugoniot  shock-energy  relation  governs,  given  by 


Eh-ph(V0-V)/2 


(2) 


then  eqns  (1)  and  (2)  may  be  combined,  eliminating  Eh,  to  yield 

p-Ehr-phl  1-(V0-V)/(2\|/)]  .  (3) 

With  the  Hugoniot  as  the  equation-of-state  reference  function,  this  form  of  the  equation  of  state 
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is  defined  only  where  the  ph  function  is  defined — namely,  between  the  reference  volume,  V0,  and 
a  limiting  compressive  volume  at  which  the  Hugoniot  asymptotes  to  infinite  pressure.  The 
thermodynamic  study  of  material  behavior  along  the  principal  Hugoniot  is,  in  and  of  itself, 
important,  sincfc  the  vast  majority  of  experimental  shock  data  obtained  lie  along  this  curve. 
Additionally,  because  most  hydrocode  implementations  of  the  equation  of  state  are  given 
explicitly  by  the  Hugoniot-based  eqn  (3),  a  pedagogical  study  of  this  equation  form  is  further 
warranted. 

We  introduce  the  nondimensional  variable  z  to  denote  the  quantity 

z~(V0-V)/( 2\|f)  ,  (4) 

which  appears  on  the  right-hand  side  of  eqn  (3).  In  terms  of  more  traditional  equation-of-state 
variables,  z  may  be  expressed  as 


z-Tp/2  ,  (5) 

where  the  compression,  p,  is  defined  as  p  -  (V0  -  V  )/V.  This  variable,  z,  will  prove  to  be  a 
convenient  term  grouping  in  which  to  express  subsequent  results.  Whereas  y\r  may  be  thought 
of  as  a  specific  volume,  scaled  by  a  factor  of  (1/T),  the  variable,  z,  on  the  other  hand,  may  be 
thought  of  as  the  nondimensional  compression,  scaled  by  a  factor  of  (172).  Thus,  in  terms  of  z, 
the  Hugoniot-based  equation  of  state  is  given  as 

p-E/y  -pA(l -z)  .  (6) 

Note  that,  under  the  assumption  that  r«r(V),  the  right-hand  side  of  this  equation  remains  a 
function  of  volume  only.  Within  the  range  of  validity  of  the  model  and  data,  eqn  (6)  may  be 
used  to  model  all  valid  thermodynamic  states  of  a  material,  on  or  away  from  the  Hugoniot 
reference  curve. 
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2.  The  Behavior  of  z 

We  may  learn  something  of  the  limiting  behavior  of  z  by  studying  eqn(6).  At  large 
enough  compressions,  the  Hugoniot  will  asymptote  vertically  toward  infinite  pressure  [since 
electronic  effects  cannot  be  modeled  within  the  F(V )  assumption].  At  the  specific  volume 
corresponding  to  this  Hugoniot  asymptote,  call  it  Vx ,  the  zero-temperature  isotherm  is  a  well- 
behaved,  finite-valued  function — thus,  pc(VJ  and  EC(VX)  are  bounded  (where  the  subscript  “c” 
denotes  conditions  along  the  zero-temperature  isotherm).  Similarly,  r*  and  thus  ip*  are  nonzero 
and  finite,  as  well.  Examining  eqn  (6)  at  the  zero  temperature  state  then,  the  left-hand  side  of 
the  equation  is  finite  and  bounded.  With  ph  at  this  specific  volume  being  infinite,  the  equation 
can  only  hold  true  if  z-»  1  at  the  specific  volume  associated  with  the  asymptotic  state  of  the 
Hugoniot.  On  the  other  hand,  at  the  reference,  unshocked  state,  where  V  equals  V0,  the 
compression  and  thus  the  value  of  z  are  both  identically  zero. 

Thus,  for  the  limits  of  unshocked  and  infinitely  shocked  states,  the  value  of  the  z  variable 
takes  on  values  of  0  and  1,  respectively.  For  shock  strengths  in  between  these  asymptotic  states, 
Segletes  [1]  showed,  with  his  Mode  II  criterion,  that  0  <  T  <  2/p  is  a  requirement  for 
thermodynamic  stability  of  any  material  for  which  T  -  T(V).  Expressed  in  terms  of  z,  this 
relation  is  simply 


0<z<l  ,  (7) 

showing  that,  indeed,  the  value  of  z  at  Vx  and  V0  bounds  the  behavior  of  z  between  Vx  and  V0. 

We  now  inquire  as  to  whether  the  variable  z  behaves  monotonically  between  these 
limiting  conditions.  We  proceed  by  differentiating,  with  respect  to  V,  the  value  of  z  given  in 
eqn  (4),  which  yields 


z'  -  -(1  A|0  (1/2  +  z  V)  •  (8) 

Since  \|r  and  z  are  always  positive  for  shock  compression,  the  behavior  of  z7  hinges  on  that  of 
\j/.  Segletes  previously  showed  [2]  that,  as  a  requirement  to  avoid  so-called  Rayleigh-slope 
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instabilities,  the  criterion  \j/y  >  — 1  needs  to  be  satisfied.  The  violation  of  this  criterion  can 
produce,  for  shocks  originating  at  elevated  temperatures,  a  postshock  pressure  below  the  preshock 
pressure— clearly,  a  nonsensical  result.  Adherence  to  the  \|/>-l  criterion  would  guarantee  the 
monotonicity  of  z  (i.e.,  z'  <  0)  for  the  domain  0  <  z  <  1/2.  Segletes  later  showed  [3]  that,  from 
the  consideration  of  isobaric  expansions  at  high  pressure,  a  stricter  condition,  "V|r/  >  0,  was 
applicable  for  r(V)  materials.  This  stricter  condition,  in  light  of  eqn(8),  would  guarantee  z- 
monotonicity  at  all  specific  volumes,  and  is  a  particularly  relevant  constraint  in  the  domain 
1/2<z<1,  which  corresponds  to  conditions  of  high  compression  and  temperature  along  the 
Hugoniot.  Thus,  under  the  assumption  of  a  material  for  which  T  -  T(V ),  to  which  we  limit  this 
discussion,  the  monotonicity  of  z  follows  from  results  previously  established  [2, 3],  Though  this 
monotonicity  is  established  for  r(V)  materials  over  the  complete  Hugoniot,  the  practical 
limitations  of  the  r(V)  assumption  (for  which  these  results  were  derived)  is  limited,  for  many 
solids,  to  a  pressure  range  of  several  megabars. 


3.  Isentropic  Behavior  Near  the  Hugoniot 

To  begin  our  examination  of  thermodynamic  behavior  along  the  Hugoniot,  we  begin  with 
the  isentrope.  Starting  with  the  shock  form  of  the  equation  of  state,  eqn  (6),  differentiate  along 
an  isentrope,  to  obtain 


(dp/dv)s  =  -phz‘  +  p'hiX-z) 


\\r{dE/dv)s  -  E\f 

\jf2 


(9) 


Knowing  that  (dE/dV  )s  -  p,  one  may  evaluate  eqn  (9)  along  the  Hugoniot  (i.e.,  along  an  isentrope 
where  it  intersects  with  the  Hugoniot),  such  that  p  and  E  take  on  values  of  ph  and  Eh, 
respectively,  while  the  expression  (dp/dV  )sh  is  used,  here,  to  denote  the  slope  of  the  isentrope 
through  a  point  on  the  Hugoniot.  One  may  eliminate  the  Hugoniot  internal-energy  term  by 
substituting  the  shock  energy  relation,  eqn  (2),  expressed  in  z  as  Eh  -  ph\\rz,  to  obtain  the  result 


(ap/dv)jA  =  p'h(\-z)  -  ^(1  +  zV  *  \ tz')  . 


(10) 
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The  product  differentiation  rule  may  be  employed  to  combine  the  terms  (z\|/  +  \|tz/)  into  (z^)7. 
The  product,  is  simply  (V0  -  V  )/2,  and  thus  its  derivative  is  -1/2.  Simplification  gives 

This  equation  relates  the  slope  of  an  isentrope  at  a  point  on  the  Hugoniot  to  the  slope  of 
the  Hugoniot  and  the  value  of  the  Hugoniot  pressure.  To  interpret  this  equation,  so  as  to  better 
understand  the  variable  z,  multiply  by  -1  and  substitute  the  definition  of  z,  eqn  (4),  into  the  last 
term  to  obtain 


-(dp/av)jA  =  - p'h(1-z. ) 


(V0~V)  z 


(12) 


The  left-hand  side  is  the  slope,  in  (p,V)  space,  of  the  isentrope  through  a  point  on  the  Hugoniot. 
It  is  expressed  as  a  linear  combination,  in  z,  of  two  terms:  the  slope  of  the  Hugoniot  itself  at  that 
point,  and  the  slope  of  the  Rayleigh  line  through  the  point  and  the  Hugoniot  origin.  The  negative 
signs  were  used  to  keep  the  physical  magnitudes  of  these  slopes  positive.  Using  the  terms  77, 
7,  and  R  to  respectively  represent  the  slopes  of  the  Hugoniot,  isentrope,  and  Rayleigh  line  through 
the  given  Hugoniot  point  (see  Figure  1),  the  expression  may  be  restated  as 

7  -  77(1  -z)  +  Rz  .  (13) 

Because  z  remains  between  0  and  1,  the  value  of  7  may  be  thought  of  as  a  weighted 
average  of  77  and  R.  Not  only  does  this  equation,  therefore,  reinforce  our  understanding  that  the 
isentropic  slope  on  the  principal  Hugoniot  must  lie  between  the  Rayleigh  and  Hugoniot  slopes 
(i.e.,  77>7>7?  >  0,  used  by  Segletes  [1]  to  prove  his  Mode  II  criterion),  but  it  also  demonstrates 
that  the  variable  z  is  intimately  related  to  these  slopes  and  to  the  disturbance  velocities  that  these 
slopes  engender.  For  example,  at  large  compressions,  where  the  Hugoniot  slope,  77,  approaches 
unbounded  values,  ideal  gas  limits  tell  us  that  the  isentropic  slope,  7,  at  a  given  volume,  can  be 
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shown  proportional  to  pressure.  The  Rayleigh  slope,  R,  also  has  this  pressure  proportionality, 
by  definition.  Thus,  in  the  high-pressure  limit,  even  though  ph  becomes  unbounded,  the  ratio  HR 
remains  finitely  bounded.  The  only  way  that  this  condition  can  hold  with  H  becoming 
unbounded  is  that,  not  only  must  z— >1,  but  z  must  functionally  approach  unity  such  that  H(l-z) 
remains  bounded,  as  well.  By  contrast,  at  the  low  end  of  the  shock-compression  spectrum 
(z«l),  eqn  (13)  indicates  that  the  isentrope  is  approximated  by  the  Hugoniot — a  well  known  fact, 
but  one  which  is  made  functionally  clear  from  the  structure  of  eqn  (13)  knowing,  further,  that  at 
z*0,  R~H. 


4.  Thermodynamic  Relations  in  Terms  of  H  and  R 

An  important  relationship  relating  the  isentropic  slope,  /,  to  the  Hugoniot  (characterized 
by  H  and  R),  in  terms  of  the  z  variable,  was  derived  in  eqn  (13).  Realizing  that  slopes  R  and  H 
are  specified  by  a  given  Hugoniot  point  and  its  (p,V)-slope,  an  intrinsic  linkage  between  the 
Hugoniot  and  Griineisen  function  becomes  undeniably  apparent.  Other  thermodynamic  paths 
through  points  on  the  Hugoniot  may  be  similarly  expressed  in  terms  of  H,  R,  and  z.  Consider 
the  thermodynamic  change  along  a  infinitesimal  path  of  arbitary  slope  -m,  in  (p,V)  space  [i.e., 
define  the  path  such  that  (dp/dV)m--nt].  It  may  be  shown  from  eqn  (6)  that 

-  W.V  *  (3B3V).  =  -pj 2  -  p,V(l  -z)  *  0 >-*>,)¥'  •  (14) 

At  points  on  the  principal  Hugoniot,  the  last  term  becomes  identically  zero.  Were  one  to  study 
material  states  off  of  the  principal  Hugoniot,  this  last  term  of  eqn  (14)  should  be  retained.  We 
restrict,  however,  the  current  analysis  to  the  Hugoniot  and  thus  discard  the  term.  Employing  the 
definitions  of  H  and  R  used  earlier,  and  substituting  -m  for  the  (p,V)-slope  in  question,  eqn  (14) 
becomes,  along  the  Hugoniot, 

my  +  {dE/dv)m  =  [H(l  -z)  -  Rz]y  .  (15) 

This  equation  is  actually  of  a  general  form,  from  which  a  result  like  eqn  (13)  may  likewise  be 
obtained.  This  equation  may  be  used  to  study  not  only  volume  derivatives  of  pressure,  but 
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volume  derivatives  of  energy,  as  well.  In  addition  to  using  it  to  study  the  Hugoniot,  isentrope, 
and  Rayleigh-line  paths,  we  will  also  use  it  to  study  the  isobaric  and  isoenergy  paths,  whose 
negative  slopes  we  denote  with  P  and  respectively.  The  following  additional  relationships  are 
readily  established  from  eqn(15)  (original  relationships  repeated  for  completeness): 


H  -  -p'  ,  (16) 

R-PkHYo-V)  ,  d7> 

I-H(l-z)  +  Rz  ,  (18) 

P-0,  (19) 

&-H(l-z)-Rz  ,  (20) 

(dE/dv)  H  =  E'h=  [-Hz-  Rz]y  ,  (21) 

(dE/dv)g  =  [if(l  ~z)  -R(l+z)]y  ,  (22) 

(dE/dv),  =  -ph  =  [~2Rz\v  ,  (23) 

{dE/dv)p  =  [H(l  -z)  -Rz]y  ,  and  (24) 

{dE/dv) f  =  0  .  (25) 


Various  graphical  interpretations  may  be  inferred  from  these  relations.  In  some  cases,  the 
interpretations  may  even  be  independent  of  z.  For  example,  eqns(21)  and  (23)  may  be 
respectively  expressed  as 


-(dE/ dv)H  =  (H  +  R)  (V0  -  V)/2  ,  and  (26) 

-(dE/dV),  =  R(V0-  V)  .  (27) 

Since  the  physical  unit  of  these  equations  is  pressure,  these  two  relations  may  be  graphically 
interpreted  in  (p,V)  space.  Figure  2  depicts  the  graphical  construction.  In  the  figure,  R' 
represents  the  line  going  through  the  Hugoniot  point  of  interest,  but  with  a  slope  negative  to  that 
of  line  R  (i.e.,  the  mirror  image  of  R  about  the  constant-pressure  line  />=/>*).  The  dE/dV 
quantities  shown  in  the  figure  are  magnitudes  only,  as  the  actual  quantities  carry  a  negative  value. 
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Fig.  1.  The  slopes  H,  7,  and  R,  shown  at  a  given  volume,  V,  define  the  respective  slopes  of 
the  Hugoniot,  isentrope,  and  Rayleigh  line  through  a  given  point  on  the  principal 
Hugoniot 


Fig.  2.  A  graphical  interpretation  of  (BE/ dV)H  compared  to  that  of  (dE/dV)„  constructed 
from  the  principal-Hugoniot  curve. 
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5.  Observations  and  Conclusions 

It  has  been  shown  that  the  state  variable,  z,  equal  to  Tp/2,  embodies  important 
thermodynamic  information.  On  the  surface,  it  represents  a  measure  of  compression,  which  has 
been  nonlinearly  scaled  by  the  factor  172.  Such  a  metric  seems  a  natural  variable  for  discussion 
of  shock  compression.  For  the  weakest  of  shocks,  its  value  approaches  zero.  In  the  strong-shock 
limit,  it  approaches  unity.  Geometrically,  it  may  be  employed  to  relate  the  properties  of 
Hugoniots  to  those  along  other  thermodynamic  paths.  Analytically,  it  behaves  under  constraints 
that  are  not  material-  or  loading-path  specific.  The  value  of  z,  coupled  with  point  and  slope 
information  on  the  Hugoniot,  provide  enough  information  to  define  the  thermodynamics  of  a 
material.  However,  the  algebraic  simplicity  by  which  thermodynamic  changes  along  paths  that 
cross  the  Hugoniot  may  be  characterized,  given  in  eqns  (18)— (25),  is  startling  and,  in  some  cases, 
amenable  to  graphical  interpretation. 

Though  the  Griineisen  assumption  [T  -  F(V)],  under  which  the  present  results  were 
derived,  can  no  longer  be  expected  to  hold  at  extremely  large  pressures  and  temperatures,  it  has 
historically  been  shown  adequate  for  many  materials  well  into  the  megabar  shock-pressure  range. 
Furthermore,  the  equation-of-state  implementation  in  virtually  all  hydrocodes  is  based  upon  the 
T(V)  assumption.  Thus,  its  study  is  of  more  than  just  academic  interest.  Because  of  the 
monotonicity  of  z,  as  shown  in  the  current  work,  required  of  materials  for  which  T  -  T(V),  a  one- 
to-one  correspondence  must  exist  between  z  and  specific  volume,  V.  Such  a  correspondence,  as 
a  minimum,  opens  up  the  possibility  that  equations  of  state  may  be  expressed  more  conveniently 
in  terms  of  z,  rather  than  V.  There  are  some  advantages  in  so  doing.  First,  since  the  Griineisen 
function  is  already  embedded  in  the  variable  z,  expressing  an  equation  of  state  in  terms  of  z 
would  effectively  recouple  the  Griineisen  function,  T,  back  to  the  equation-of-state  reference 
curve  (e.g.,  the  Hugoniot  or  cold-curve).  Such  a  coupling  is  important,  as  it  has  been  shown  [1] 
that  arbitrary  selection  of  a  Griineisen  functional  form  for  a  given  reference  Hugoniot  can  often 
lead  to  thermodynamic  inconsistency.  By  establishing  a  functional  coupling  between  the 
reference  curve  and  the  Griineisen  function  (even  if  it  were  empirically  derived),  one  could  check 
the  functional  form  for  thermodynamic  consistency,  rather  than  exhaustively  testing  every 
empirical  combination  of  reference  curve  and  Griineisen  function. 
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Of  possible  greater  significance  is  the  fact  that  the  introduction  of  the  variable  z  creates, 
perhaps,  a  framework  for  the  development  of  a  generic  equation  of  state  for  solids  in  terms  of 
primitive,  macroscopic  variables.  This  framework  is  possible  because  not  only  is  z  unit- 
dimensionless,  but  for  any  arbitrary  material  subject  to  shock  loading,  the  value  of  z  will  proceed 
from  zero  to  unity  as  shock  strength  increases.  There  are  many  approaches  by  which  an  equation 
of  state  in  z  might  be  formulated.  One  could  specify  the  reference  curve  in  terms  of  z  directly, 
as  in  ph(z)  or  Eh(z).  Such  direct  specification  though  would  likely  be  the  result  of  an  empirical 
fitting  procedure.  One  might  alternately  attempt  to  specify  the  functional  coupling  by 
interrelating  slopes  H  and  R  as  a  function  of  z.  This  latter  approach  is  functionally  equivalent 
to  relating  ph'  to  ph  by  way  of  z,  or  alternately  Eh'  to  Eh  by  way  of  z.  Such  techniques  are  more 
amenable  to  analytical  or  semi-analytical  derivation,  since  there  do  exist  constraints  which  govern 
thermodynamic  functions  and  their  derivatives.  Of  course,  the  challenge  here  is  in  finding  the 
appropriate  constraint  that  precisely  maps  the  derived  equation  of  state  onto  the  experimentally 
observed  data. 
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SOCORRO  NM  87801 

3  AT  .1 1  ANT  TECHSYSTEMS  INC 
ATTN  T  HOLMQUIST  MN11  2720 
RSTRYK 

G  R  JOHNSON  MN11  2925 
600  SECOND  ST  NE 
HOPKINS  MN  55343 


1  MARVIN  L  ALME 

2180  LOMA  LINDA  DR 
LOS  ALAMOS  NM  87544-2769 

1  APPLIED  RESEARCH  ASSOC  INC 
ATTN  JEROME  D  YATTEAU 

5941  S  MIDDLEFIELD  RD  SUITE  100 
LITTLETON  CO  80123 

2  APPLIED  RESEARCH  ASSOC  INC 
ATTN  DENNIS  GRADY 
FRANK  MAESTAS 

4300  SAN  MATEO  BLVD  SE 
ALBUQUERQUE  NM  87110 

1  BATTELLE 

ATIN  ROBER  M  DUGAS 

7501  S  MEMORIAL  PKWY  SUITE  101 

HUNTSVILLE  AL  35802-2258 

3  BOEING  AEROSPACE  CO 
SHOCK  PHYSICS  &  APPLIED  MATH 
ENGINEERING  TECHNOLOGY 
ATIN  R  HELZER 

T  MURRAY 
J SHRADER 
PO  BOX  3999 
SEATTLE  WA  98124 

1  BOEING  HOUSTON  SPACE  STN 
ATTN  RUSSELL  F  GRAVES 
BOX  58747 
HOUSTON  TX  77258 

1  BRIGS  CO 

ATIN  JOSEPH  E  BACKOFEN 
2668  PETERSBOROUGH  ST 
HERNDON  VA  20171-2443 

1  CALIFORNIA  RSCH  &  TECHNOLOGY 
ATIN  M  MAJERUS 
PO  BOX  2229 
PRINCETON  NJ  08543 

1  CENTURY  DYNAMICS  INC 
ATTN  N  BIRNBAUM 
2333  SAN  RAMON  VALLEY  BLVD 
SAN  RAMON  CA  94583-1613 
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1  COMPUTATIONAL  MECHANICS 
CONSULTANTS 
ATTN  J  A  ZUKAS 
PO  BOX  11314 
BALTIMORE  MD  21239-0314 

1  CYPRESS  INTERNATIONAL 
ATTN  A  CAPONECCHI 
1201  E  ABINGDON  DR 
ALEXANDRIA  VA  22314 

1  DEFENSE  TECHNOLOGY  INTL.  INC 
ATTN  D  E  AYER 
THE  STARK  HOUSE 
22  CONCORD  ST 
NASHUA  NH  03060 

1  DESKIN  RESEARCH  GROUP  INC 
ATTN  EDWARD  COLLINS 
2270  AGNEW  RD 
SANTA  CLARA  CA  95054 

3  DOW  CHEMICAL  INC 
ORDNANCE  SYSTEMS 
ATTN  CHANEY 
A  HART 
B  RAFANIELLO 
800  BUILDING 
MIDLAND  MI  48667 

1  G  E  DUVALL 

5814  NE  82ND  COURT 
VANCOUVER  WA  98662-5944 

3  DYNA  EAST  CORP 
ATTN  P  C  CHOU 
R  CICCARELLI 
WFLIS 

3620  HORIZON  DRIVE 
KING  OF  PRUSSIA  PA  19406 

3  DYNASEN 

ATTN  JACQUES  CHAREST 
MICHAEL  CHAREST 
MARTIN  LILLY 
20  ARNOLD  PL 
GOLETA  CA  93117 


1  R  J  EICHELBERGER 
409  W  CATHERINE  ST 
BEL  AIR  MD  21014-3613 

1  ELORET  INSTITUTE 

ATTN  DAVID  W  BOGDANOFF  MS  230  2 
NASA  AMES  RESEARCH  CENTER 
MOFFETT  FIELD  CA  94035 

3  ENIG  ASSOCIATES  INC 
ATTN  J  ENIG 
D  J  PASTINE 
M  COWPERTHWATIE 
SUITE  500 

11120  NEW  HAMPSHIRE  AVE 
SILVER  SPRING  MD  20904-2633 

1  EXPLOSIVE  TECHNOLOGY 
ATTN  M  L  KNAEBEL 
PO  BOX  KK 
FAIRFIELD  CA  94533 

1  GB  TECH  LOCKHEED 
ATTN  JAY  LAUGHMAN 
2200  SPACE  PARK  SUITE  400 
HOUSTON  TX  77258 

2  GB  TECH  LOCKHEED 

ATTN  LUCILLE  BORREGO  C23C 
JOE  FALCON  JR  C23C 
2400  NASA  ROAD  1 
HOUSTON  TX  77058 

6  GDLS 

38500  MOUND  RD 
ATTN  W  BURKE  MZ436-21-24 
G  CAMPBELL  MZ436-30-44 
D  DEBUSSCHER  MZ436-20-29 
J  ERIDON  MZ436-21-24 
W  HERMAN  MZ  435-01-24 
S  PENTESCU  MZ436-21-24 
STERLING  HTS  MI  48310-3200 

2  GENERAL  RESEARCH  CORP 
ATTN  A  CHARTERS 
TMENNA 
PO  BOX  6770 

SANTA  BARBARA  CA  93160-6770 
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2  GRC  INTERNATIONAL 

ATTN  TIMOTHY  M  CUNNINGHAM 
WILLIAM  M  ISBELL 
5383  HOLLISTER  AVE 
SANTA  BARBARA  CA  93111 

6  INST  OF  ADVANCED  TECHNOLOGY 
UNIVERSITY  OF  TX  AUSTIN 
ATTN  STEPHEN  J  BLESS 
JAMES  CAZAMIAS 
HARRY  D  FAIR 
THOMAS  M  KIEHNE 
DAVID  LITTLEFIELD 
MIKE  NORMANDIA 
4030-2  W  BRAKER  LN 
AUSTIN  TX  78759 

1  INTERNATIONAL  RESEARCH  ASSOC 
ATTN  D  L  ORPHAL 
4450  BLACK  AVE 
PLEASANTON  CA  94566 

1  INTERPLAY 

ATTN  F  E  WALKER 
18  SHADOW  OAK  RD 
DANVILLE  CA  94526 

1  R  JAMESON 
624  ROWE  DR 
ABERDEEN  MD  21001 

1  KAMAN  SCIENCES  CORP 
ATTN  DENNIS  L  JONES 
2560  HUNTINGTON  AVE  SUITE  200 
ALEXANDRIA  VA  22303 

8  KAMAN  SCIENCES  CORP 
ATTN  J  ELDER 
RICHARD  P  HENDERSON 
DAVID  A  PYLES 
FRANK  R  SAVAGE 
JAMES  A  SUMMERS 
JAMES  S  WILBECK 
TIMOTHY  W  MOORE 
THY  YEM 

600  BLVD  S  SUITE  208 
HUNTSVILLE  AL  35802 


3  KAMAN  SCIENCES  CORP 
ATTN  SHELDON  JONES 
GARY  L  PADEREWSKI 
ROBERT  G  PONZINI 

1500  GRDN  OF  THE  GODS  RD 
COLORADO  SPRINGS  CO  80907 

4  KAMAN  SCIENCES  CORP 
ATTN  NASIT  ARI 
STEVE  R  DIEHL 
WILLIAM  DOANE 
VERNON  M  SMITH 

PO  BOX  7463 

COLORADO  SPRINGS  CO  80933-7463 

1  D  R  KENNEDY  &  ASSOC  INC 
ATTN  D  KENNEDY 
PO  BOX  4003 

MOUNTAIN  VIEW  CA  94040 

1  KERLEY  PUBLISHING  SERVICES 
ATTN  G  I  KERLEY 

PO  BOX  13835 

ALBUQUERQUE  NM  87192-3835 

2  KTECH  CORPORATION 
ATTN  FRANK  W  DAVIES 
LARRY  M  LEE 

901  PENNSYLVANIA  NE 
ALBUQUERQUE  NM  87110 

1  LIVERMORE  SOFTWARE  TECH  CORP 
ATTN  J  O  HALLQUIST 
2876  WAVERLY  WAY 
LIVERMORE  CA  94550 

1  LOCKHEED  MARTIN  MISSLE  &  SPACE 
ATTN  WILLIAM  R  EBERLE 
PO  BOX  070017 
HUNTSVILLE  AL  35807 

3  LOCKHEED  MARTIN  MISSILE  &  SPACE 
ATTN  M  A  LEVIN  ORG  81  06  BLDG  598 
M  R  MCHENRY 

T  A  NGO  ORG  81  10  BLDG  157 
111  LOCKHEED  WAY 
SUNNYVALE  CA  94088 
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4  LOCKHEED  MISSILE  &  SPACE  CO 
ATTN  JOHN  R  ANDERSON 
WILLIAM  C  KNUDSON 
S  KUSUMI  0  81  11  BLDG  157 
J  PHILLIPS  0  54  50 
PO  BOX  3504 
SUNNYVALE  CA  94088 

1  LOCKHEED  MISSILE  &  SPACE  CO 
ATTN  R  HOFFMAN 
SANTA  CRUZ  FACILITY 
EMPIRE  GRADE  RD 
SANTA  CRUZ  CA  95060 

1  LOCKHEED  NASA  JSC 
SPACE  SCIENCE  BRANCH 
ATTN  JAMES  HYDE 
BOX  58561  MC  B22 
HOUSTON  TX  77258 

1  LOCKHEED  MARTIN  AEROSPACE 
ATTN  D  R  BRAGG 
PO  BOX  5837  MP  109 
ORLANDO  FL  32855 

1  MCDONNELL  DOUGLAS 
ASTRONAUTICS  CO 
ATTN  B  L  COOPER 
5301  BOLSA  AVE 
HUNTINGTON  BEACH  CA  92647 

2  ORLANDO  TECHNOLOGY  INC 
ATTN  DANIEL  A  MATUSKA 
MICHAEL  GUNGER 

PO  BOX  855 

SHALIMAR  FL  32579-0855 

1  PHYSICAL  SCIENCES  INC 
ATTN  PETER  NEBOLSINE 
20  NEW  ENGLAND  BUS  CTR 
ANDOVER  MA  01810 

3  PHYSICS  INTERNATIONAL 
ATTN  R  FUNSTON 

G  FRAZIER 
L  GARNETT 
PO  BOX  5010 
SAN  LEANDRO  CA  94577 


NO.  OF 

COPIES  ORGANIZATION 

1  PRC  INC 

ATTN  J  ADAMS 

5166  POTOMAC  DR  #103 

KING  GEORGE  VA  22485-5824 

1  RAYTHEON  ELECTRONIC  SYSTEMS 
ATTN  R  KARPP 
50  APPLE  HELL  DRIVE 
TEWKSBURY  MA  01876 

1  ROCKWELL  INTERNATIONAL 
ROCKETDYNE  DIVISION 
ATTN  H  LEIFER 
16557  PARK  LN  CIRCLE 
LOS  ANGELES  CA  90049 

1  ROCKWELL  MISSILE  SYS  DIV 
ATTN  T  NEUHART 
1800  SATELLITE  BLVD 
DULUTH  GA  30136 

1  SAIC 

ATTN  MICHAEL  W  MCKAY 
10260  CAMPUS  POINT  DR 
SAN  DIEGO  CA  92121 

1  SHOCK  TRANSIENTS  INC 
ATTN  DAVID  DAVISON 
BOX  5357 

HOPKINS  MN  55343 

2  SIMULATION  &  ENG  CO  INC 
ATTN  ELSA  I  MULLINS 
STEVEN  E  MULLINS 

8840  HWY  20  SUITE  200  N 
MADISON  AL  35758 

2  SOUTHERN  RESEARCH  INSTITUTE 
ATTN  LINDSEY  A  DECKARD 
DONALD  P  SEGERS 
PO  BOX  55305 

BIRMINGHAM  AL  35255-5305 


NO.  OF  NO.  OF 

COPIES  ORGANIZATION  COPIES  ORGANIZATION 


5  SRI  INTERNATIONAL 
ATTN  JAMES  D  COLTON 
D CURRAN 
RKLOOP 
RL  SEAMAN 
D  A  SHOCKEY 
333  RAVENSWOOD  AVE 
MENLO  PARK  CA  94025 

2  TELEDYNE  BROWN  ENGR 
ATTN  JIM  W  BOOTH 
MARTIN  B  RICHARDSON 
PO  BOX  070007  MS  50 
HUNTSVILLE  AL  35807-7007 

1  ZERNOW  TECHNICAL  SVCS  INC 
ATTN  LOUIS  ZERNOW 
425  W  BONITA  AVE  SUITE  208 
SAN  DIMAS  CA  91773 


ABERDEEN  PROVING  GROUND 

20  DIR,  USARL 

ATTN:  AMSRL-WM,  I  MAY 
AMSRL-WM-MC,  J  WELLS 
AMSRL-WM-PA,  S  HOWARD 
AMSRL-WM-PB,  A  ZIELINSKI 
AMSRL-WM-PC,  R  PESCE-RODRIGUEZ 
AMSRL-WM-PD,  G  GAZONAS 
AMSRL-WM-T,  W  F  MORRISON 
AMSRL-WM-TA, 

M  BURKINS 
W  GILLICH 
W  BRUCHEY 
JDEHN 
GFILBEY 
W  A  GOOCH 
HW  MEYER 
E  J  RAPACKI 
J  RUNYEON 
AMSRL-WM-TB, 

R  FREY 
P BAKER 
R  LOTTERO 
J  STARKENBERG 


30  DIR,  USARL 

ATTN:  AMSRL-WM-TC, 

W  S  DE  ROSSET 
T  WBJERKE 
R  COATES 
F  GRACE 
KKIMSEY 
M  LAMPSON 
D  SCHEFFLER 
S  SCHRAML 
GSILSBY 
B  SORENSEN 
R  SUMMERS 
W  WALTERS 
AMSRL-WM-TD, 

SCHOU 
A  M  DIETRICH 
D  DANDEKAR 
K  FRANK 
J  HARRISON 
M  RAFTENBERG 
ARAJENDRAN 
M  SCHEIDLER 
S  SCHOENFELD 
S  SEGLETES  (5  CP) 
J  WALTER 
T  WRIGHT 
AMSRL-WM-WD, 

J  POWELL 
APRAKASH 
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1 


1 


1 


1 


1 


1 


ORGANIZATION 

AERONAUTICAL  &  MARITIME 
RESEARCH  LABORATORY 
ATTN  N  BURMAN 
S  CIMPOERU 
D  PAUL 
PO  BOX  4331 
MELBOURNE  VIC  3001 
AUSTRALIA 

EMBASSY  OF  AUSTRALIA 
ATTN  R  WOODWARD 
COUNSELLOR  DEFENCE  SCIENCE 
1601  MASSACHUSETTS  AVE  NW 
WASHINGTON  DC  20036-2273 

ABTEILUNG  FUER  PHYSKALISCHE 
CHEMIE 

MONTANUNIVERSITAET 
ATTN  E  KOENIGSBERGER 
A  8700  LEOBEN 
AUSTRIA 

PRB  S  A 

ATTN  M  VANSNICK 

AVENUE  DE  TERVUEREN  168  BTE  7 

BRUSSELS  B  1150 

BELGIUM 

ROYAL  MILITARY  ACADEMY 
ATTN  E  CELENS 
RENAISSANCE  AVE  30 
B1040  BRUSSELS 
BELGIUM 

BULGARIAN  ACADEMY  OF  SCIENCES 
SPACE  RESEARCH  INSTITUTE 
ATTN  VALENTIN  GOSPODINOV 
1000  SOFIA  PO  BOX  799 
BULGARIA 

CANADIAN  ARSENALS  LTD 
ATTN  P  PELLETIER 
5  MONTEE  DES  ARSENAUX 
VILLIE  DE  GRADEUR  PQ  J5Z2 
CANADA 

DEFENCE  RSCH  ESTAB  SUFFIELD 
ATTN  D  MACKAY 

RALSTON  ALBERTA  TOJ  2N0  RALSTON 
CANADA 


NO.  OF 

COPIES  ORGANIZATION 

1  DEFENCE  RSCH  ESTAB  SUFFIELD 
ATTN  CHRIS  WEICKERT 
BOX  4000  MEDICINE  HAT 
ALBERTA  TIA  8K6 
CANADA 

1  DEFENCE  RSCH  ESTAB  VALCARTIER 
ARMAMENTS  DIVISION 
ATTN  R  DELAGRAVE 
2459  PIE  XI  BLVD  N 
PO  BOX  8800 

CORCELETTE  QUEBEC  GOA  1R0 
CANADA 

1  UNIVERSITY  OF  GUELPH 
PHYSICS  DEPT 
ATTN  C  G  GRAY 
GUELPH  ONTARIO 
NIG  2W1 
CANADA 

1  CEA 

ATTN  ROGER  CHERET 
CEDEX  15 

313  33  RUE  DE  LA  FEDERATION 

PARIS  75752 

FRANCE 

1  CEA 

CISI  BRANCH 
ATTN  PATRICK  DAVID 
CENTRE  DE  SACLAY  BP  28 
GIF  SUR  YVETTE  91192 
FRANCE 

1  CEA/CESTA 

ATTN  ALAIN  GEILLE 
BOX  2  LE  BARP  33114 
FRANCE 

6  CENTRE  D’ETUDES  DE  GRAMAT 
ATTN  SOLVE  GERARD 
CHRISTIAN  LOUPIAS 
PASCALE  OUTREBON 
J  CAGNOUX 
C  GALLIC 
J  TRANCHET 
GRAMAT  46500 
FRANCE 
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2  CENTRE  D’ETUDES  DE  LMEIL-VALENTON 
ATIN  CHRISTIAN  AUSSOURD 
JEAN-CLAUDE  BOZIER 

SAINT  GEORGES  CEDEX 
VILLENEUVE  94195 
FRANCE 

3  CENTRE  D’ETUDES  DE  VAUJOURS 
ATTN  PLOTARD  JEAN-PAUL 
ERIC  BOTTET 

TAT  SIHN  VONG 
BOriE  POSTALE  NO  7 
COUNTRY  77181 
FRANCE 

6  CENTRE  DE  RECHERCHES 
ET  D’ETUDES  D’ARCUEEL 
ATTN  D  BOUVART 
C  COTTENNOT 
S  JONNEAUX 
HORSINI 
S  SERROR 
FTARDIVAL 

16  BIS  AVENUE  PRIEUR  DE 
LA  COTE  D’OR 
F94114  ARCUEEL  CEDEX 
FRANCE 

1  DAT  ETBS  CETAM 

ATTN  CLAUDE  ALTMAYER 
ROUTE  DE  GUERRY  BOURGES 
18015 
FRANCE 

1  ETBS  DSTT 

ATIN  P  BARNIER 
ROUTE  DE  GUERAY 
BOITE  POSTALE  712 
18015  BOURGES  CEDEX 
FRANCE 

1  FRENCH  GERMAN  RESEARCH  INST 
ATIN  CHANTERET  P-Y 
CEDEX  12  RUE  DE  I’INDUSTRIE 
BP  301 

F68301  SAINT-LOUIS 
FRANCE 


5  FRENCH  GERMAN  RESEARCH  INST 
ATIN  HANS-JURGEN  ERNST 
FRANCIS  JAMET 
PASCALE  LEHMANN 
KHOOG 
HLERR 

CEDEX  5  5  RUE  DU  GENERAL 
CASSAGNOU 
SAINT  LOUIS  68301 
FRANCE 

1  LABORATOIRE  DE  TECHNOLOGIE  DES 
SURFACES 

ECOLE  CENTRALE  DE  LYON 
ATTN  VINETP 
BP  163 

69131  ECULLY  CEDEX 
FRANCE 

1  BATTELLE  ENGENIEUTECHNIK  GMBH 
ATIN  W  FUCHE 
DUESSELDORFFER  STR  9 
ESCHBORN  D  65760 
GERMANY 

1  CONDAT 

ATIN  J  KIERMEIR 
MAXMILIANSTR  28 
8069  SCHEYERN  FERNHAG 
GERMANY 

1  DEUTSCHE  AEROSPACE  AG 
ATIN  MANFRED  HELD 
POSTFACH  13  40 
D  86523  SCHROBENHAUSEN 
GERMANY 

1  DIEHL  GBMH  AND  CO 
ATIN  M  SCHILDKNECHT 
HSCHBACHSTRASSE  16 
D  90552  ROTBENBACH  AD  PEGNITZ 
GERMANY 
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5  ERNST  MACH  INSTITUT 
ATTN  VOLKER  HOHLER 
E  SCHMOLINSKE 
E  SCHNEIDER 
ASTILP 
KTHOMA 
ECKERSTRASSE  4 
D-7800  FREIBURG  I  BR  791  4 
GERMANY 

1  EUROPEAN  SPACE  AGENCY 
ATTN  WALTER  FLURY 
ESOC5 

ROBT  BOSCHT  STRASSE 
DARMSTADT  6100 
GERMANY 

3  FRAUNHOFER  INSTITUT  FUER 
KURZZEITDYNAMIK 
ERNST  MACH  INSTITUT 
ATTN  H  ROTHENHAEUSLER 
H  SENF 

E  STRASSBURGER 
HAUPTSTRASSE  18 
D79576  WEIL  AM  RHEIN 
GERMANY 

3  FRENCH  GERMAN  RESEARCH  INST 
ATTN  HARTMUTH  F  LEHR 
ROLF  HUNKLER 
ERICH  WOLLMANN 
POSTFACH  1260 
WEIL  AM  RHEIN  D-79574 
GERMANY 

2  IABG 

ATTN  M  BORRMANN 
H  G  DORSCH 
EINSTEINSTRASSE  20 
D  8012  OTTOBRUN  B  MUENCHEN 
GERMANY 

1  INGENIEURBURO  DEISENROTH 
AUF  DE  HARDT  33  35 
D5204  LOHMAR  1 
GERMANY 


3  TU  CHEMNITZ-ZWICKAU 
ATTN  I  FABER 
L  KRUEGER 
LOTHAR  MEYER 

FAKULTAET  FUER  MASCHINENBAU  U. 

VERFAHRENSTECHNIK 
SCHEFFELSTRASSE  110 
09120  CHEMNITZ 
GERMANY 

1  TECHNISCHE  UNIVERSITAT  MUENCHEN 
ATTN  EDUARD  B  IGENBERGS 
RICHARD  WAGNER  STR  18  111 
MUENCHEN  2  D8000 
GERMANY 

1  BHABHA  ATOMIC  RESEARCH  CENTRE 
HIGH  PRESSURE  PHYSICS  DIVISION 
ATTN  N  SURESH 
TROMBAY  BOMBAY  400  085 
INDIA 

5  RAFAEL  BALLISTICS  CENTER 
ATTN  EREZ  DEKEL 
YEHUDA  PARTOM 
G  ROSENBERG 
Z  ROSENBERG 
Y  YESHURUN 
PO  BOX  2250 
HAIFA  31021 
ISRAEL 

1  TECHNION  INST  OF  TECH 
FACULTY  OF  MECH  ENGNG 
ATTN  SOL  BODNER 
TECHNION  CITY 
HAIFA  32000 
ISRAEL 

1  IHI  RESEARCH  INSTITUTE 
STRUCTURE  &  STRENGTH 
ATTN:  TADASHI  SHIBUE 
1-15,  TOYOSU  3 
KOTO,  TOKYO  135 
JAPAN 
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1  ESTEC  CS 

ATTN  DOUGLAS  CASWELL 
BOX  200  NOORDWUK 
2200  AG 
NETHERLANDS 

2  EUROPEAN  SPACE  AGENCY  ESTEC 
ATTN  LUCY  BERTHOUD 
MICHEL  LAMBERT 

POSTBUS  BOX  299  NOORDWUK 

NL2200  AG 

NETHERLANDS 

4  PRINS  MAURITS  LABORATORY 
ATTN  H  J  REITSMA 
EDWARD  VAN  RIET 
HPASMAN 
R  YSSELSTEIN 
TNO  BOX  45 
RUSWIJK  2280AA 
NETHERLANDS 

1  ROYAL  NETHERLANDS  ARMY 
ATTN  J  HOENEVELD 

V  D  BURCHLAAN  31 
PO  BOX  90822 

2509  LS  THE  HAGUE 
NETHERLANDS 

4  HIGH  ENERGY  DENSITY  RESEARCH  CTR 
ATTN  VLADIMIR  E  FORTOV 
GENADH  I  KANEL 

V  A  SKVORTSOV 
O  YU  VOJOBIEV 
IZHORSKAJA  STR  13/19 
MOSCOW  127412 
RUSSIAN  REPUBLIC 

1  INSTITUTE  OF  CHEMICAL  PHYSICS 
ATTN  A  YU  DOLGOBORODOV 
KOSYGIN  ST  4  V  334 
MOSCOW 

RUSSIAN  REPUBLIC 


3  INSTITUTE  OF  CHEMICAL  PHYSICS 
RUSSIAN  ACADEMY  OF  SCIENCES 
ATTN  A  M  MOLODETS 
S  V  RAZORENOV 
A  V  UTKIN 

142432  CHERNOGOLOVKA 
MOSCOW  REGION 
RUSSIAN  REPUBLIC 

3  INSTITUTE  OF  MECH  ENGINEERING  PROBLEMS 

ATTN  V  BULATOV 
D  INDEITSEV 
Y  MESCHERYAKOV 
BOLSHOY,  61,  V.O. 

ST  PETERSBURG  199178 
RUSSIAN  REPUBLIC 

2  IOFFE  PHYSICO  TECHNICAL  INSTITUTE 
DENSE  PLASMA  DYNAMICS 
LABORATORY 

ATTN  EDWARD  M  DROBYSHEVSKI 
A  KOZHUSHKO 
ST  PETERSBURG  194021 
RUSSIAN  REPUBLIC 

1  IPERAS 

ATTN  A  A  BOGOMAZ 
DVORTSOVAIA  NAB  18 
ST  PETERSBURG 
RUSSIAN  REPUBLIC 

2  LAVRENTYEV  INST.  HYDRODYNAMICS 
ATTN  LEV  A  MERZHIEVSKY 
VICTOR  V  SILVESTROV 
NOVOSIBIRSK  630090 

RUSSIAN  REPUBLIC 

1  MOSCOW  INST  OF  PHYSICS  &  TECH 
ATTN  S  V  UTYUZHNIKOV 
DEPT  OF  COMPUTATIONAL 
MATHEMATICS 
DOLGOPRUDNY  1471700 
RUSSIAN  REPUBLIC 


1  RESEARCH  INSTITUTE  OF  MECHANICS 
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